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Topological meaning of Z 2 numbers in time reversal invariant systems 
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We show that the Z2 invariant, which classifies the topological properties of time reversal 
invariant insulators, has deep relationship with the global anomaly. Although the second Chern 
number is the basic topological invariant characterizing time reversal systems, we show that 
the relative phase between the Kramers doublet reduces the topological quantum number Z to 
7, 2 . 

KEYWORDS: Quantum spin Hall effect, Time reversal symmetry, Z2 topological invariant, global anomaly 



As one of possibilities for spintronics, spin Hall effect 
(SHE) has been proposed, 1,2 and much effort has been 
made to realize it experimentally. 3-5 In particular, the 
quantum spin Hall effect (QSHE) has been attracting 
much current interest due to its unique topological char- 
acteristics. 6-8 Although it may be regarded as a spin ver- 
sion of integer quantum Hall effect (IQHE), their differ- 
ences with respect to topological properties have been 
revealed. It is well known that IQHE states are classified 
by integers which are topological invariant called first 
Chern number. 9-12 Kane and Mele 6 proposed that time 
reversal invariant states are labeled by Z 2 numbers, and 
much work has been done to clarify its meaning. 13-20 The 
Z 2 number implies that time reversal invariant systems 
have two phases: One is the ordinary insulator, whereas 
the other is the topological insulator 21 ' 22 in which QSHE 
is expected. 

On the other hand, two decades ago, Avron et. al. 23 
argued that the basic topological invariant for time re- 
versal invariant systems is the second Chern number. Qi, 
Hughes and Zhang recently proposed a generic topologi- 
cal field theory in 4+1 dimensions, where basic topolog- 
ical properties are classified by the second Chern num- 
ber. 24,25 From this, they have constructed topological 
insulators in two and three spatial dimensions by the use 
of the dimensional reduction. 

In this paper, we address the question what is the fun- 
damental relationship between the second Chern number 
and Z2 number. We show that time reversal invariant sys- 
tems are indeed labeled by the second Chern numbers. 
However, taking account of the relative phase ambiguity 
between a Kramers doublet, the topological classification 
by Z falls into two sectors labeled by Z 2 . This reflects 
7r 4 (Sp(l))=Z 2 (more generically, 7r 4 (Sp(n))=Z 2 ), and it 
has close relationship with the global anomalies found by 
Witten. 26,27 

In generic Hamiltonian with broken time reversal sym- 
metry, the eigenvalue degeneracies have codimension 
three. 23 This means that due to level crossings, Hamilto- 
nian depending on three parameters could have some iso- 
lated "monopoles" in the parameter space, and on some 
two dimensional subspaces (typically the 2-sphere, S 2 ), 
the first Chern number takes a nontrivial integer if the 
monopoles exist inside the sphere. On the other hand, 



in the time reversal invariant systems with T 2 = — 1, 
the eigenvalue degeneracies have codimension five. 23,25 
Therefore, to reveal the topological properties of such 
systems, we need to examine Hamiltonians with five pa- 
rameters. The minimal model may be the quaternionic 
2x2 hermitian matrix or equivalently, 4x4 hermitian ma- 
trix with time reversal invariancc, which indeed depends 
on five parameters. 23,25 

Let H(x) be the Hamiltonian satisfying TH{x)T^ 1 = 
H(x), where x denotes a set of five parameters x a (a = 
1, ••• ,5). Each eigenstate is doubly-degenerate, which 
is referred to as Kramers doublet. Such an eigenstate 
with the minimum twofold degeneracy may be called 
quaternionically simple. 23 At certain points in five di- 
mensional (5D) parameter space, level crossings among 
several Kramers doublets occur. Then, we can consider 
a specific 4D subspace in which all the eigenstates are 
quaternionically simple, that is, there are no degeneracies 
other than the Kramers degeneracy. We will show that 
topological numbers can be defined on this 4D space, 
which we assume to be S 4 below. 25 

Owing to the assumption above, Kramers doublets are 
distinguishable on S 4 , and therefore, we can examine 
each Kramers doublet separately. Let ^ be a set of eigen- 
functions of a specific Kramers doublet: 25 



One finds the time reversal constraint, 
T* = * J, 

where the matrix J is defined by 

1 



J 



-1 



(1) 



(2) 



(3) 



which operates on the space of the Kramers doublet 
Berry's gauge potential is defined by 



.4 



(4) 



where d is the exterior derivative with respect to x a , 
d = dx a d/dx a . The time reversal constraint (2) ensures 
that A = —J~ 1 A T J , implying that A € sp(l) algebra. 
In the case of generic 2n Kramers multiplet, which may 
happen if the Hamiltonian has some other symmetries, 
the gauge potential belongs to sp(n). Therefore, the first 
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Chcrn number vanishes since tr F = holds identically, 
where F is the field strength 2-form defined by F = 
dA+A 2 . Arvon et. al. 23 showed that for the time reversal 
invariant systems, the basic topological invariant is the 
second Chern number, 

In the case of the minimum quaternionic 2x2 model men- 
tioned above, they derived C2 = ±1 for two Kramers 
doublets. 

The nontrivial Chern number is due to the obstruc- 
tion of the gauge fixing. 10, 25, 28 Suppose that we need two 
patches to span smooth eigenfunctions over S . Typical 
example is such that S 4 is divided into two hemispheres 
which are topologically equivalent to the 4D discs, de- 
noted here as D± . In the overlap region n D A _ , the 
two kinds of the wave functions \P±, defined on D±, re- 
spectively, are related with each other through the gauge 
transformation (the transition function): 

= *_fif. (6) 

The time reversal constraint (2) imposes the condition 
g = J~ x g*J, which tells that g eSp(l). The gauge po- 
tentials A± , defined on the patches D± , respectively, are 
thus related with each other as 

A + =g- 1 A_g + g- 1 dg. (7) 

As in the case of the first Chcrn number, the sec- 
ond Chcrn number c 2 can be written by the transi- 
tion function g in the following manner: First, note that 
tr F 2 = duj 3 , where ojs(A,F) is the Chern-Simons 3- 
form uj 3 (A,F) = tr (AF - A 3 /3). It is not gauge in- 
variant; its transformation law is given by lu 3 (A 9 , F 9 ) = 
u 3 (A,F) + da 2 - {l/3)tr (g^dg) 3 , where A 9 and F 9 de- 
note the gauge transforms of A and F by a generic g, 
and ai2 is a certain 2-form which vanishes after the inte- 
gration. Thus, we find that the second Chern number C2 
defined in eq. (5) is given by 

C2 ^/ s3 tr(rt) 3 , (8) 

where S 3 denotes 3-sphere as the boundary of D 4 ,. This 
tells that C2 is given by the winding number of the tran- 
sition function g in eqs. (6) and (7) over S 3 . 

So far we have discussed that the Chern number (5) 
would be one of natural topological numbers for the spec- 
ification of time reversal invariant systems. 23 In more 
generic cases with 2n Kramers multiplct, g belongs to 
Sp(n), and therefore, the Chern number is also given by 
a winding number of the transition function g G Sp(n) on 
S 3 . This is ensured by 7r 3 (Sp(n))=Z (n > 1). All these 
suggest that time reversal invariant states would be clas- 
sified by integers like the IQHE. However, it contradicts 
the Z2 classification of time reversal invariant insulators 
by Kane and Mele, 6 who assert that there are only two 
phases, trivial insulators and topological insulators in 
two spatial dimensions. Their argument is quite clear, es- 
pecially from the point of view of the edge states, which 
naturally pick up the topologically nontrivial property 
of the bulk: 11 Even number of edge states are unstable 
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against time reversal invariant perturbations. Then, how 
does the Chern number classification, which corresponds 
to the bulk picture, match the Z2 classification? To an- 
swer this question, we should note that in the above cal- 
culations we have used the wave function of the Kramers 
doublet defined in eq. (1). This equation means that once 
we choose a wave function tp for one state, we fix that 
of the other paired state as —Tip. However, we can basi- 
cally attach an arbitrary phase to the paired state. This 
reminds us of the pfaffian formula proposed by Fu and 
Kane, 15 in which the relative phase between the Kramers 
doublet plays an important role. Therefore, it is very cru- 
cial to take into account of this relative phase ambiguity. 
To be more precise, this phase is described as follows: A 
transformation ip — > e %e l 2 i\) induces 

r{6) = e"^/ 2 )- 3 . (9) 

For fixed 9, this never affect the topological properties 
described above. Indeed, r(9) is a gauge transformation 
belonging to Sp(l), under which the Chern number (5) 
is invariant. However, if the phase 9 is changed adiabati- 
cally and moved along a loop, it reveals a new topological 
property of the transition function which is different from 
the second Chern number. As we will see below, such a 
topological aspect of the transition function is nonpertur- 
bative. Note that r(0) = 1 and r(2ir) = —1. The overall 
sign of ivP has nothing to do with the relative phase, so 
that the period of 9 should be < 9 < 2tt. Eq. (9) yields 
the transformation for the transition function, 

g(x) - r- 1 (9)g(x)r(9) = g r (x, 9). (10) 

The r.h.s above denotes the "rotated" transition function 
which now depends on 9 as well as x. Note that the period 
of g r with respect to 9 is indeed 2ir, g r (x, 0) = g r {x, 2ir). 
Thus, the transition function can be regarded as a func- 
tion defined on the 4D space S 3 x S 1 , where the latter S 1 
is the loop spanned by 9. Under the 2ir rotation induced 
by the adiabatic change of 9, the transition function falls 
into two classes, as demonstrated by Witten. 26 ' 27 

Without use of algebraic topology, Witten 26,27 in- 
vented a method of calculating the global anomaly in 
terms of the perturbative anomaly. This method is very 
helpful, also in the present problem, to define and to 
compute the Z2 numbers. In what follows, by the use of 
the technique developed by Witten, we first define the 
Z2 number and next show that it is related to the second 
Chcrn number C2 such that odd (even) C2 is nontriv- 
ial (trivial) element of Z 2 . To this end, we embed the 
Sp(l)~SU(2) transition function g into SU(3) such that 

9( X )=( 9{X) (ID 

Likewise, r{9) is embedded into SU(3) as 
f{9) = e -*(e/2)(A 3 -v^A 8 ) 

= diag(l,e* e- i9 ), (12) 

where A a is the standard Gell-Mann matrices. Here, we 
have used the fact that the special form g in eq. (11) is 
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invariant under a particular SU(3) transformation e raA 
with A 8 = diag(l, 1,— 2)/V3 for any real a. Although a 
natural embedding may be only the A 3 = diag(l, —1,0) 
part in eq. (12), this invariance enables us to introduce 
another rotation with respect to A 8 , as in eq. (12). Fur- 
thermore, we introduce a radial parameter p (0 < p < 1) 
into f to enlarge the parameter space S 3 x S 1 to five 
dimensional disc- like space S 3 x D 2 , whose boundary 
(p = 1) is S 3 x S\ i.e, d(S 3 x D 2 ) = S 3 x S 1 . At the 
boundary we assume f(9, 1) to be eq. (12), whereas in- 
side S 3 x D 2 (0 < p < 1), f(9 1 p) is a generic element of 
SU(3). The following is a convenient extrapolation from 
the boundary to p = proposed by Witten; 



1 



f(9,p) 



pe 



pe~ i9 



(13) 



Note that r(9, p = 0) is independent of 9, implying that 
the space spanned by (8, p) has indeed a disc geometry. 
Thus, the rotated transition function g r in (10) is em- 
bedded into SU(3) as 

g r (x,9,p)=f- 1 (9,p)g(x)r(9,p). (14) 

Now the Z 2 number D can be defined by 



D = T[g r ], 



where 



(15) 



(16) 



240tt 2 J S 3 XD 2 

To be precise, D in eq. (15) should be defined mod 2tt, as 
will be clarified momentarily. In what follows, we study 
the properties of D thus defined. 

First, we note that it is topological, since T[g r ] is in- 
variant under infinitesimal changes of g r , T[g r + Sg r ] = 
T[g r ], provided that g r satisfies the boundary condition 
(11) and (12). We next show that the Z 2 number D, in 
the case of the embedding above, is related to the second 
Chern number. To this end, we note the relation 



r[gh} = T[g]+T[h} + AT[g,h}, 



(17) 



for generic SU(3) fields g and h. After the integration 
over p, the last term is given by 



AT[g,h] = 



48tt 2 



tr 



d(S 3 xD 2 ) 



{g^dgfdhh- 1 



(18) 



+ g^dgidhh- 1 ) 3 + -(g^dgdhh- 1 ) 2 

Applying the identity (17) to eq. (16), we obtain T[g r ] = 
AT [f -1 ,<7f] + Ar[<7,f]. By making use of the fact that 



r(9, p) in (13) is independent of 9 at p 
perform the integration over 9, 



0, we can further 



trdff- 1 [{dg-r 1 ? + {9- l d~g)^ 



* JS 3 xS^ 



= 7TC 2 . (19) 

It thus turns out that D defined by eq. (15) is propor- 
tional to the second Chern number. In more generic cases 



with 2n Kramers multiplct, the above relation (19) still 
holds by embedding Sp(n) into SU(3n). 

Finally, we study the modulo 2n property of D. Basi- 
cally, it is due to the dependence of the way of embedding 
in the present formulation. To show this, it is convenient 
to introduce another function which denotes the relative 
change of g r from g, 

G(x,9)=g r (x,9)g- 1 (x). (20) 

This function helps us to see a more precise origin of 
the Z 2 number. As we have noted, x a and 9 span S 3 
and S 1 , respectively, and therefore, g r is a mapping from 
S 3 x 5* 1 to Sp(l). On the other hand, the new mapping 
(20) satisfies the relation, G(x,0) = G{x,2ir) = 1. This 
boundary condition as well as the fact that 7Ti(Sp(l))=0 
make it possible to regard S 3 x S 1 as S 4 , and thus, the 
mapping G should be classified by 7T4(Sp(l))=Z 2 . In a 
similar way above, the mapping (20) can be embedded 
into SU(3) as 

G{x,8,p) =g r (x,9,p)g- 1 (x,0,p), (21) 

where g r has been defined in eq. (14). The boundary 
condition is G(x,0,p) = G(x,2ir,p) — 1. Therefore, we 
can define the Z 2 number as 



D = T[G], mod 2tt, 



(22) 



where T is the same as in eq. (16) but with the integra- 
tion domain D 5 whose boundary is S 4 mentioned-above, 
dD 5 = S 4 . The modulo 2ir property is now clear: If D 5 
and D 5 are two different discs but with the same bound- 
ary S 4 , then we have 

T D .[G\-T D , 6 [G] 



■ I {G^dGf. 



240** ' ' (23) 

The r.h.s. is manifestly a multiple of 2ir, which is en- 
sured by tt 5 (SU(3))=z! It is known 26 - 27 - 29 that when G 
is a trivial element of 7T4(Sp(l)), its embedding G gives 
D — 0, whereas in the case with nontrivial G it yields 
D = it. Now, we will show below that two D's defined in 
eqs. (15) and (22) indeed coincide. This leads to the con- 
clusion that g(x) on S 3 is classified by Z (second Chern 
number), but even (odd) elements can be continuously 
deformed into different even (odd) elements through the 
relative phase 9 between the Kramers doublet. The non- 
trivial element in 7T4(Sp(l)), which should describe the 
nontrivial topological states with time reversal symme- 
try, is thus given by g(x) with odd second Chern num- 
bers. We believe that this is the origin of Z 2 character- 
istics of the topological number specifying the phases of 
time reversal invariant systems. 

At first sight, the two definitions (15) and (22) seem 
different, but straightforward calculations lead us to 
T[G] = T[g r }. To show this, we note that eq. (17) yields 

T[G] = T[g r ] + AT, (24) 

where the second term is essentially given by eq. (18) 
with g — > g r (x, 9, p) and h — > g^ 1 (x, 0, p), 



AT : 



48tt 2 



L 



tr 



S 3 xS 1 



+ 9r ld 9r(g 1 dg) 



(9r ld 9rfg ^ 



-Ag r 1 dg r g 1 dg) 2 



(25) 
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In what follows, we show that this term vanishes, Ar = 
0, and hence, T[G] — T[g r ] holds. To this end, we first 
note that the I-form g~ x dg is written as 



9 1 dg 







(26) 



where v is sp(I)-valucd I-form v = g~ 1 dg. Then, one 
finds that the rotated 1-form can be decomposed into 
two parts: 



g r 1 dg r =u + v r 



(27) 

The former is the I-form u including only d9, while the 
latter is rotated v with dx only, 

u = f~ 1 g~ 1 rdr~ 1 gr + r~ 1 dr 
— 1 \3zz 



g- L ygr-\ 3 )d6 



ud6 



f 1 vf 



v r 



(28) 



Kramers doublet has an ambiguity in relative phase be- 
tween the wave functions of the doublet. Through the 
adiabatic change of this phase, the wave function with a 
certain Chern number can be continuously deformed into 
other wavefunctions with different Chern numbers. What 
is important is that in this deformation the evenness and 
the oddness of the Chern numbers are conserved. We be- 
lieve that this is the origin of the Z 2 property of the 
topological number for the time reversal invariant sys- 
tems. Mathematically, it can be described by the fourth 
homotopy class of the transition function with the rela- 
tive phase degree of freedom, 7r 4 (Sp(n))=Z2. 
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1 {9 1t3 9 - r3 ) r = u ° ra an d v r 



where u — |r~ 
Now eqs. (26) and (28) tell that the fields in the integrand 
in eq. (25) can be written only by sp(l)-valued fields, 

r 2ir 



Ar = 



48tt 2 



dS 



/ tru(v 2 v 
Js 3 



+ vv 2 + v 3 — v r vv r — vv r v) 



(29) 

The first three terms in the integrand vanish under the 
trace because sp(l) is pseudo-real. The last two terms 
seem to be complicated, because 9 and x are coupled 
together. Nevertheless, 9 is basically decoupled from x, 
and therefore, we can carry out the integration over 9. 
To see this, let us denote v — v a r a 1 and calculate, for 
example, the following term; 

v r vv r = r~ 1 vrvr~ 1 vr 

= v 1 v 2 v 3 e abc r~ 1 T a rT b r~ 1 T c r 



v 1 v 2 v 3 t 1 t 2 t 3 2(r 2 + r~ 2 + 1) 
|,V+r- 2 + l). 



(30) 



i d9 [ tTu[v 3 (r 2 +r- 2 + l)] . (31) 

y Jo Js 3 



The same formula holds for the other term in eq. (29), 
and we finally arrive at 

r 2ir 

i 

A1 = 48tt2 3! Jo j s3 

Since v 3 is proportional to the identity matrix, we find 
that the terms with r 2 and r~ 2 vanish by integration 
over 9. The same is true for the last term because of the 
identity trr a w 3 = 0. Thus, we conclude that Ar = 0, 
and therefore, T[G] = T[g r ]. It follows from eq. (19) that 
for trivial and nontrivial mapping G — r~ 1 grg~ 1 , g has 
even and odd second Chern number C2, respectively. 

In summary, we have shown that the second Chern 
number is relevant topological number for the time re- 
versal invariant systems with T 2 = —1. However, the 
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